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We consider a two dimensional Kondo lattice model with 
exchange J and hopping t in which three out of four impurity 
spins are removed in a regular way. At the particle-hole sym- 
metric point the model may be studied with auxiliary field 
quantum Monte Carlo methods without sign problems. To 
achieve the relevant energy scales on finite clusters, we in- 
troduce a simple method to reduce size effects by up to an 
order of magnitude in temperature. In this model, a metal- 
lic phase survives up to arbitrarily low temperatures before 
being disrupted by magnetic fiuctuations which open a gap 
in the charge sector. We study the formation of the heavy- 
electron state with emphasis on a crossover scale T* defined 
by the maximum in the resistivity versus temperature curve. 
The behavior of thermodynamic properties such as specific 
heat as well as spin and charge uniform susceptibilities are 
studied as the temperature varies in a wide range across T* . 
Within our accuracy T* compares well to the Kondo scale of 
the related single impurity problem. Finally our QMC resuls 
are compared with mean-field approximations. 
PACS numbers: 71.27.-fa, 71.10.-w, 71.10.Fd 



I. INTRODUCTION 

Kondo insulators as well as heavy fermion materials 
[Q,^ are believed to be described by Kondo lattice mod- 
els virith antiferromagnetic exchange J and hopping t. 
In this framework, the Kondo insulator - exemplified by 
Ce3Bi4Pt3 - corresponds to a special band filling in which 
there is exactly one conduction electron per impurity 
spin. The origin of the charge and spin gaps is easy to un- 
derstand in the limit of large J/ 1 where each conduction 
electron is trapped by an impurity spin to form a Kondo 
singlet. Starting from this point there are two ways of 
generating metallic states which will described the heavy 
electron state. On one hand, one can remove conduction 
electrons thus leaving uncompensated impurity spins. In 
the limit J/t —>■ oo the problem maps onto the U —>■ oo 
Hubbard model where the Kondo singlets are represented 
by empty sites and the bachelor - or unscreened - spins by 
electrons [1). Thus, a metallic state is expected. On the 
other hand one can keep the number of conduction elec- 
trons constant and deplete the lattice of impurity spins. 
Starting from Ce3Bi4Pt3 this amounts to replacing Ce by 
La Each missing impurity spin liberates a conduc- 
tion electron. As can be seen easily in the strong coupling 
limit this conduction electron is localized since its motion 
involves the breaking of adjacent Kondo singlets. Hence 



a bound state with magnetic properties appears within 
the gap Q. If the depletion is large enough the induced 
localized levels may overlap coherently to form a metallic 
state. Alternatively this metallic state where there are 
more conduction electrons than impurity spins may be 
viewed as evolving from dilute impurity limit as realized 
for example in Cea;Lai — xCug Note however that in 
the small (large) x limit the Ce (La) atoms are randomly 
substituted for La (Ce) in the LaCug (CeCug) crystal 
structure. 

In the next section, we first consider within a mean- 
field approximation both above describe routes to 
obtain metallic states. The paramagnetic solution of 
this mean-field approximation is nothing but the large-N 
(N corresponds to the SU(N) symmetry of the impurity 
spins) saddle point ||^-[lO||. Here, we concentrate on a 
special choice of depletion scheme for the localized spins. 
Three out of four impurity spins are removed in a regular 
way thus ensuring that in the limit J/t —^ oo the ground 
state is metallic. At the mean-field level both routes to 
obtain a metallic state (doping and above describe deple- 
tion scheme) lead to the same metallic ground state: a 
Fermi-liquid with arbitrarily large effective mass (or ar- 
bitrarily low coherence temperature) as the coupling J/ 1 
is reduced. The mean-field solution has several pitfalls. 
For instance, the imaginary part of the self-energy van- 
ishes at all temperatures, the Kondo scale corresponds 
to a phase transition rather than to a crossover scale and 
charge fluctuations on impurity sites are allowed since 
constraints are imposed only on average. 

For depleted Kondo lattices it is possible to carry out 
auxiliary field quantum Monte Carlo (QMC) simulations 
without encountering the infamous sign problem. This 
comes with the restriction of particle-hole symmetry for 
the conduction band. Nevertheless this gives us the op- 
portunity to compare exact results with mean-field ap- 
proximations. There are two challenges to circumvent 
before achieving good results. The first one is the imple- 
mentation of the exchange interaction J along with the 
constraint of singly occupied impurity sites. This prob- 
lem has been solved in pd[ | and extensively discussed in 
Jl2t . The second challenge lies in controlling finite size 
effects stemming from the conduction electrons. This is 
extremely important in the framework of Kondo physics 
since screening takes place on an energy scale set by the 
Kondo temperature Tk which is exponentially small in 
the weak coupling limit. If the level spacing of the con- 
duction electrons exceeds this small energy sca le, fin ite 
size effect will dominate the result. In section [II B we 
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present a simple method to reduce by up to an order mag- 
nitude in temperature the onset of size effects in thermo- 
dynamic quantities (specific heat, spin and charge sus- 
ceptibiUtics) for a given lattice size. 

Section [V summarizes our numerical results. Ground 



state, thermodynamic and transport properties are com- 
puted on lattice sizes up to 500 sites. We show that our 
simulations reproduce the typical form of the resistiv- 
ity versus temperature curve observed in heavy fermion 
materials: an initial rise as a function of decreasing 
temperature originating from the onset of screening of 
the magnetic impurities then followed by a sharp down- 
turn signaling the onset of coherence between individual 
screening clouds. Energy scales as a function of coupling 
strength are computed. The last section is devoted to 
conclusions and discussion of our results. 



II. MODELS AND MEAN-FIELD 
APPROXIMATION 

Our starting point is the two-dimensional Kondo lat- 
tice model (KLM): 

Here i runs over the lattice sites of a two dimensional 
square lattice: i = nSx + may, n,m : 1 ■ ■ ■ L with Sx, dy 
the lattice vectors. The set of lattice points {R} denote 
the positions of impurity spins, cl creates a conduction 

of 



electron on site i with z-component of spin a and S — 

R 

\ Y.a,a' f\ ^^'T.'^'/fl.cr' ^^"^ ^^'^ ^P^^^ ^/^ Operators with a 
the Pauli spin matrices. A similar definition holds for 
5^. Since the KLM stems from the strong coupling limit 
of the periodic Anderson model ||l^ the charge degrees 
of freedom on the f-sites are frozen. Hence, a constraint 
of one electron per /-site holds for KLM. 

In the notation of Eq. |l| the usual KLM model is ob- 
tained when R = na^ + may (see Fig. ^). In principle, 
the depletion scheme of impurity spins can be random - 
a situation which would reflect the experimental realiza- 
tion - or regular. As a first step we have chosen a regu- 
lar depletion and considered R — n [2ax] + m[ax + 2ay] 
(see Fig. ^). From now onwards we will refer to this 
model as the depleted Kondo lattice model (DKLM) . The 
unit cell of the DKLM has five orbitals: one localized for 
the impurity spin and four delocalized for the conduction 
electrons. The Bravais lattice is spanned by the lattice 
vectors: ai = 2dx and a2 = dx + 2a„. 
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FIG. 1. (a) The standard KLM. The square lattice is 
spanned by the lattice vectors a^, Sy. (b) The depleted KLM 
(DKLM). The corresponding Bravais lattice is spanned by 
2ax and + 2%. Each unit cell contains a single localized 
orbital for the impurity spin and four delocalized orbitals ac- 
commodating the conduction electrons. In both (a) and (b) 
the top layer denotes the conduction electrons with hopping 
matrix element t. The bottom layer are the magnetic impuri- 
ties which couple to the conduction electrons via the exchange 
J. (c) the Brillouin zone of the DKLM. The reciprocal lattice 
vectors read bi = ^(0, 1) and 62 = ^(1, -1/2) 

Having defined the model, we now address the ques- 
tion: to which degree does the KLM a,t ric < 1 have 
the same physics as the DKLM at = 4? Here ric 
corresponds to the number of conduction electrons per 
unit cell. As mentioned previously, in the strong cou- 
pling limit the KLM with band filling ric maps onto the 
U ^ 00 Hubbard model with band filling he — I — ric [|| . 
The nature of the metallic state is not clear due to the 
strong correlations. However for fillings far from the dop- 
ing induced Mott metal-insulator transition, ric = 1, it 
is reasonable to believe that beyond the one-dimensional 
case the ground state is a Fermi liquid as described for 
example by a Gutzwiller approximation. Due to the spe- 
cial choice of the depletion, the strong coupling DKLM 
at ric = 4 is a Fermi liquid. In this limit, each impurity 
spin will capture a conduction electron and the remnant 
conduction electrons will form the metallic state. To be 
more precise at J/t 00 we obtain a three band model, 
where a single conduction band crosses the Fermi level, 
6^ = 0. The thus obtained strong coupling Fermi line is 
shown in Fig. |^c. 

In the weak coupling limit, we can address the above 
question within a mean-field approximation . One can 
write the KLM and DKLM models as: 



m{R})^j2<H^%. 
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where the last term accounts precisely for the spin 

flip terms: ^ Y^r (/i.^/fla^a'^^.T + 4,i^«.T4,T^«.i) ' 
With the above rewriting of the KLM we can introduce 
order parameters both for magnetism and Kondo screen- 
ing. 



-vire^^'^ and 



1 / t _ t \ 

(ft cs , + cl / s ,) = (ft eg , + ct fs , ) = -V. (3) 

where Q correspond to the relevant wave vector for mag- 
netic order. Concentrating in a first step on paramag- 
netic solutions TTis = m/ = we obtain the mean- field 
Hamiltonian 
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where corresponds to the chemical potential, A to a 
Lagrange multiplier enforcing the constraint of single 
occupancy on the /-sites on average and Nimp is the 
number of impurity spins. The saddle point equations, 
^ = ^ = 0, "/T'f^ = where F is the free energy 
may then be solved to obtain the mean-field solution. 
The thus obtained saddle point corresponds precisely to 
that obtained in the large-N approximation [p|-JlO[. 
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FIG. 2. Paramagnetic mean-field solution of the KLM at 
ric = 0.5 (a) The band structure. Energies are measured with 
respect to the Fermi energy. The inset corresponds to the 
quasiparticle weight of the conduction band, (b) Fermi line, 
(c) Optical conductivity. Note that there is a Drude peak 
at a; = 0. (d) Drude weight normalized by the sum rule as 
well as the coherence temperature as a function of J/t. The 
inset corresponds to the specific heat from which we obtain 
the coherence temperature, (see text) 

Fig. H summarizes the paramagnetic mean-field solu- 
tion for the KLM at ric = 0.5. The band structure plotted 
in Fig. ^ consists of two bands (hybridized bands). As 
apparent, the band at the Fermi energy is flat and has 
a small quasiparticle weight (Fig. inset). The cor- 
responding Fermi surface (Fig. ^b), covers 75% of the 
Brillouin zone, so that both /—electrons and conduction 
electrons are included in the Luttinger volume. Given the 
band structure, the optical conductivity follows (Fig. ^). 
The large mass (i.e. combined small quasiparticle weight 
and flat dispersion relation ) of the conduction electrons 
at the Fermi energy leads to a small Drude weight. The 
feature at higher frequencies are nothing but interband 
transitions. The Drude weight normalized by the sum 
rule as a function of J/t is plotted in Fig. |^. At small 
values of J/t it becomes arbitrarily small reflecting the 
large mass of the charge carries. Along with the small 
Drude weight one expects a small coherence tempera- 
ture. The coherence temperature ( Tcoh ) - the energy 
scale at which the Fermi liquid character of the ground 
state becomes apparent - may be defined for example by 
the energy scale below which the specific heat is linear 
in T (inset Fig. ||d). As expected the J/t dependence 
of Tcoh tracks the Drude weight. In the mean- field ap- 
proach the Kondo temperature, Tk, corresponds to the 
temperature scale at which the order parameter V van- 
ishes. Below Tk the magnetic impurities are screened. 
In the mean-field approximation, this screening does not 
take place via the formation of a many-body state of the 
impurity spin and conduction electrons but via charge 
fiuctuations on the impurity sites. In this approxima- 
tion, charge fluctuations are allowed since the constraint 
of single occupancy of /-sites is only imposed on average. 
As apparent by the jump in the specific heat (see inset of 
Fig. ^) Tk corresponds to a a phase transition rather 
than to a crossover scale. For a recent discussion of the 
Kondo and coherence scales within this approximation, 
the reader is referred to |10|. 
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FIG. 3. Paramagnetic mean-field solution for the DKLM 
at ric — 4:. ( i.e. the unit cell contains four conduction elec- 
trons) (a) Optical conductivity, (b) Drude weight normalized 
by the sum rule, (c) Specific heat versus T/Tk- 

We now turn our attention to the DKLM. Figs. ||a,b 
plot the optical conductivity as well as the Drude weight 
versus J/t within the paramagnetic mean-field approxi- 
mation. As apparent, there is a clear similarity with the 
above presented results for the KLM. From the specific 
heat data, (Fig. || c) one can extract the coherence tem- 
perature. Here, some care has to be taken since the model 
has a van-Hove singularity at the Fermi level. Hence the 
coherence temperature corresponds to the energy scale 
below which the specific heat follows a Tln(l/r) law. 
Hence Tcoh ^ O.ITr: for the considered couphngs. 

Next, we consider magnetic degrees of freedom for the 
DKLM. As mentioned above, at J/t = oo the model 
maps onto free electrons. The underlying particle-hole 
symmetry of the model directly leads to nesting so that 
the J/t = oo point will be unstable to magnetic order- 
ing with wave vector bi/2 (Fig. |l|c). At weak couplings 
the Ruderman-Kittel-Kasuya-Yossida (RKKY) interac- 
tion favors magnetic order rather that Kondo screening 
iQ. Hence, for all values of J/t we anticipate a magnet- 
ically ordered ground state. Fig. ^ plots the mean-field 
transition temperatures corresponding to Kondo screen- 
ing, Tk, (below Tk, V ^ 0) and magnetism, Ts {ms ^ 
nif ^ 0). As apparent, magnetic order is always present. 
For values of J/t > 0.8 the mean- field Kondo scale ex- 
ceeds the magnetic scale T5. For this parameter range 
and as a function of decreasing temperature, one expects 
to be able to study the formation of the heavy electron 
state prior to the onset of magnetic fluctuations which 
due to nesting will drive the system to an insulator. The 
QMC simulations presented in section IV are aimed at 
studying this crossover. 
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FIG. 4. (a) Mean-field energy scales for the DKLM at 
Uc — 4. Tk corresponds to the Kondo scale below which 
V 7^ 0. Ts corresponds to a magnetic energy scale be- 
low which magnetic ordering with wave vector Q = fei/2 
- schematically shown in Fig. (b) - occurs. Tk and Ts 
mark phase transitions. For J/t > 0.8 both Kondo screen- 
ing {V / 0) and magnetic ordering {nis 7^ 0, m/ 7^ 0) are 
present in the ground state. At J/t < 0.8 the Kondo scale 
vanishes since in the mean-field approximation magnetic or- 
dering freezes the impurity spins. The coherence temperature 
(not plotted in the figure) lies roughly an order of magnitude 
lower that the Kondo scale. Hence for most of the considered 
couplings it will be masked by the magnetic instability. 



III. NUMERICAL TECHNIQUES 

In this section we will first briefly summarize the nu- 
merical techniques used. We then turn our attention on 
a simple method to dramatically reduce the size effects 
for the simple case of free electrons. This is the key for 
producing size independent results down to low temper- 
atures for the DKLM. However, this method is general 
and may prove useful for other applications. 



A. Quantum Monte Carlo algorithms 

We have used ground state and finite temperature ver- 
sions of the auxiliary field QMC method to investigate 
the DKLM. The T = method is based on the equation: 



(^olOl^o) (*T|e 
= iim — 

et^oo 



(*o|*o> 



KLM |\I/jn) 



(5) 



where |^t) is required to be non-orthogonal to the 
ground state |5'o)- Finite temperature properties are cal- 
culated in the grand canonical ensemble via: 



(0) 



(6) 



Details of how to implement the QMC method without 
generating a sign problem in the particle-hole filled case 
have been introduced an extensively described for the 
KLM in Refs ||ll|,|l^ . Since precisely the same tricks may 
be used for the DKLM, we refer the reader to those arti- 
cles for further details. We now concentrate on another 
important issue, size effects, which become extremely se- 
vere when one wishes to study metallic states down to 
low temperatures. 

Finally, to determine the single impurity Kondo scale 
for our given band-structure, we have used the Hirsch- 
Fye impurity algorithm [T5| ]. The same formulation as 
the lattice problem is readily implemented in this ap- 
proach. 
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B. Size effects and magnetic fields 

Reducing size effects in numerical simulations is im- 
portant, since many models of correlated electron sys- 
tems may be solved numerically only on small lattices. 
Size effects become particularly severe when the ground 
state turns out be be a metallic state with large coherence 
temperature. On the other hand, insulators are charac- 
terized by the localization of the wave function and are 
hence rather insensitive to boundary conditions on finite 
sized systems. It thus becomes apparent, that the worst 
case scenario for severe size effects are just free electrons 
in a tight binding approximation: 



H 



-t ^ clcj + H.c. 



(7) 



In many cases before turning on the interaction which 
will automatically restrict the size of the lattice under 
consideration it is important to control size effects for 
this simple case. We will concentrate on the two di- 
mensional case on a torus geometry which for the above 
model reduces to imposing periodic boundary conditions: 
ct ^ = cl, ct ^ — where L is the linear length of 

the lattice lying in the ex, Sy plane. 

To reduce size effects on thermodynamic quantities one 
may in principle consider the Hamiltonian: 



ill) 



H.c. 



(8) 



where t:^^{L) are arbitrary hopping parameters which 
have to satisfy 



hm tj-(i) 



-t. 



(9) 



Clearly this choice of hopping matrix elements on finite 
lattices will break the lattice symmetry. This is a price 
we are willing to pay provided that the convergence as a 
function of system size of thermodynamics quantities is 
greatly improved. Eq. ^ nevertheless guarantees that in 
the thermodynamic limit this symmetry is restored. To 
determine the hopping matrix elements ijj(i) so as to 
reduce size effects on say the specific heat, Cv{L,T) = 

dE(L) 
dT ' 



one may minimize 



X 



^[a(i,T)-a(i = oo,T)]' 



(10) 



The magnetic field is introduced via the Peirls phase fac- 
tors: 



H{L) 
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(11) 



with Bl{x) = V X Al{x) and $0 the flux quanta. The 
torus geometry imposes restrictions on the field. 
Since, a translation in the argument of the vector po- 
tential may be absorbed in a gauge transformation: 



Al{x + Lex)^ Al{x) + Vxx{x), 
Al{x + Ley) = Al{x) + VXy{x), 

we chose, the boundary condition 



cl 



g*"X^(i)^t gt ^ 



i+Le 

to satisfy the requirement: 

[H{L),TLsA^[H{L),TLSy] = G 



(12) 



(13) 



(14) 



Here, Tg corresponds to a translation by x. However, 
magnetic translation operators belong to the magnetic 
algebra [|l6|: 



_ {L>L^ xLey)-B 



(15) 



Thus, to obtain a single valued wave function the 
condition of flux quantization has to be satisfied: 

n where n is an integer. Here, we consider 



{Le'^xLey)-B 
5o 



a static magnetic field running along the z-axis perpen- 
dicular to the X, y plane in which the lattice lies. Hence, 
the smallest magnetic field which we can consider on a 
given lattice size satisfies: 



$0 



1. 



(16) 



With this choice of magnetic field and associated vector 
potential Eq. holds. 



where the sum extends over a given range of tempera- 
tures. Taking into account only amplitude modulations 
of the hopping matrix elements leads already to a cum- 
bersome minimization problem which does not provide 
satisfactory results. 

Instead of carrying out a complicated minimization 
problem we can try to guess which matrix elements 
tr^iL) will minimize size effects. It turns out that in- 
troducing a magnetic field produces remarkable results. 



5 



sl7$o = 1 



L = 32 













3^ 




T — I — I — I — I — I — r 

-4 -3 -2 -1 1 2 3 4 -3 -2 -1 1 2 3 4 
Uj/t Uj/t 

FIG. 5. Density of states N{lj) = jfJ2r^(^^^) 
(right column) and without (left column) magnetic 
Here, we consider 5 — O.Olt. 
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field. 



To illustrate the reduction of size effects caused by the 
inclusion of the magnetic field, we first consider the sin- 
gle particle density of states. In a basis where H{L) is 
diagonal, H{L) = E^=i en7i7« with 4 = Em7jnt^m,n 
and U^U — I, the local density of states reads: 



N{r,uj) 



N 



iS 



(17) 



where (5 is a positive infinitesimal and N the total num- 
ber of sites. Since the magnetic field breaks translation 
invariance (it is the site dependent vector potential which 
enters the Hamiltonian) N{r,uj) is site dependent. Aver- 
aging over sites yields the density of states 7V(w) plotted 
in Fig. ^. As apparent, without the magnetic field and 
up to L = 32, N{uj) is dominated by size effects for the 
considered value of 6 = O.Oli. In contrast, the presence of 
the magnetic field provides remarkable improvements. In 
particular the van-Hove singularity is well reproduced al- 
ready on L = 16 lattices and at L = 32 the result is next 
to exact for the considered value of S. It is instructive to 
look at the L = 8 case with and without magnetic fields. 
When B is turned on, the degeneracy of levels is lifted. 
Each level - apart from the e„ = level which is two fold 
degenerate - is nondegenerate. This is precisely what one 
expects for Landau levels which have degeneracy L^B/^q 
which is unity in our case. This provides an intuitive un- 
derstanding of why the method works so well. Since each 



level becomes singly degenerate, the single particle states 
cover homogeneously the the energy range of the band- 
width. Clearly this can only be achieved by breaking the 
lattice symmetry on finite sized systems. 

We now turn our attention to the specific heat coeffi- 
cient 7 ~ Cy/T (see Fig. ^a,b). As apparent, for a given 
system size, the inclusion of the magnetic field provides a 
gain of more than one order of magnitude in the temper- 
ature scale at which size effects set in. In particular the 
ln(l/T) behavior of 7 due the the van-Hove singularity 
becomes apparent already on L = 6 lattices. 



Biy^o = 1, {n) 




FIG. 6. Specific and spin susceptibility versus temperature 
without (a,c) and with (b,d) magnetic field. The curves from 
right to left correspond to increasingly large lattices as de- 
noted in the figure. The dashed line corresponds to the exact 
result. 

Upon inspection a similar improvement is obtained for 
the spin susceptibility (see Fig. ||c,d). Note that since 
we are dealing with free electrons the charge and spin 
susceptibilities are identical. 

One crucial question is whether the magnetic field will 
introduce a sign problem in the numerical simulations. It 
turns out that for DKLM at the particle-hole symmetric 
point it does not ]T^ . Hence we can incorporate it in 
our simulations and benefit from the above demonstrated 
drastic improvement in reduction of size effects. 

Other schemes have been proposed to reduce size ef- 
fects. In particular, averaging over boundary conditions 
has been suggested ||l^,|l^. This method has the advan- 
tage of not breaking translation symmetry. However, the 
averaging requires several simulations and is hence com- 
putationally expensive. In contrast with the presented 
method the improvement in reduction of size effects is 
obtained within a single simulation. 
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IV. NUMERICAL RESULTS 

Our numerical results for the DKLM are summarized 
in this section. We first concentrate on ground-state 
properties and check for magnetic ordering as a func- 
tion of J /t. Various thermodynamic quantities (specific 
heat, spin and charge susceptibilities) as well as the opti- 
cal conductivity and hence resistivity are then analyzed 
as a function of temperature at different values of J /t. 



A. Spin-spin correlations at T= 

To detect long-range magnetic ordering, we compute 
the spin-spin correlation function at T = 0: 



S{Q) = E 



(18) 



Here we consider Q — bi/2 and the sum runs over the 
impurity spins. In the presence of long-range magnetic 
order, S{Q) is an extensive quantity so that the staggered 

moment, = limL^oo \J S{Q)/L^, is finite. denotes 
the number of unit-cells. 




0.2 , 0.3 

FIG. 7. S{Q) at T = for the DKLM as a function of 
system size and various values of J/t. The L = 10 point 
corresponds to a lattice of 10 x 10 unit cells and hence 500 
orbitals. After extrapolation to the thermodynamic limit, we 
obtain the staggered moment plotted in the inset. 

Figure g plots S{Q) for various values of J/t. The 
inset in Fig. |^ plots the staggered moment obtained after 
extrapolation of S{Q) to the thermodynamic limit. The 
data is thus consistent with long-range magnetic ordering 
at least up to J/t = 2. Clearly as J/t grows and the 
staggered moment becomes small, it is increasingly hard 
to distinguish it from zero. Nevertheless, and based on 



the mean- field result, we expect the staggered moment 
to vanish only in the limit J/t — > oo. 

Following arguments put forward in Ref. ||T^ Kondo 
screening and magnetic ordering to coexist: impurities 
are partially screened and the remnant magnetic mo- 
ment orders. Depending upon the relative energy gains of 
both mechanisms a staggered moment of arbitrary mag- 
nitude can be produced. In the weak coupling limit, 
the energy gain for magnetism is the RKKY interaction. 
This interaction leads to an effective Heisenberg model 
between impurity spins with exchange Jrkky{<i) = 
— J^Rex ((f, uj — Q) where x (9j ^) corresponds to the spin 
susceptibility of the conduction electrons. On the other 
hand screening occurs on an energy scale set by the 
Kondo temperature: Tk ~ We~^^'^ where W is the 
band width. Hence for small values of J/t the system 
will gain energy by ordering magnetically. Although the 
wave vector of magnetic ordering will depend on the de- 
tails of the spin susceptibility, this source of energy gain 
is present irrespective of the band structure. On the 
other hand, the strong coupling result is an artifact of 
the model. In this limit, if Kondo screening were com- 
plete, a nested Fermi surface would be produced as in the 
J/t — » oo limit. Hence, in the ground state, the system 
will always gain energy by keeping Kondo screening par- 
tial. The remnant magnetic moment will order with the 
nesting wave vector thereby opening a gap and gaining 
energy. 

Were it not for nesting a finite critical value of the cou- 
pling J/t at which a quantum phase transition between 
ordered and disordered magnetic states would occur. Un- 
fortunately, this is beyond the scope of the QMC since 
the absence of particle-hole symmetry and hence nesting 
leads to a sign problem. 



B. Transport and thermodynamic properties. 

Here, we consider finite temperature properties of the 
model. As discussed above, the ground state is a mag- 
netic insulator irrespective of J/t. The salient features of 
the model however, lie in its thermodynamic properties. 

We first consider the real part of the optical conduc- 
tivity as obtained from the Kubo formula, (7{lu,T). This 
quantity is related to the imaginary time current-current 
correlation functions via: 



(J(t)J(O)) = J dLoK{uj,T)a{uj,T) 
with K[uj,T) 



TT 1 - e-f^'^ ' 



(19) 



Here J is the current operator along the x or y lattice 
direction and () represents an average over the finite- 
temperature ensemble. The above inverse Laplace trans- 
form, to obtain the optical conductivity is carried out 
with the ME ||l9| method. We use the paramagnetic 
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mean-field solution with substantial broadening as de- 
fault model. Note that the order parameter which 
is proportional to the spin-spin correlation between the 
conduction electron and impurity spin is estimated from 
the Monte Carlo run and not by solving the saddle point 
equations. 



0.1 0.2 ' 0.3 0.4 0.5 




FIG. 8. Optical conductivity and related resistivity as 
a function of temperature at J/t — 1.6. The data in Figs. 
(b),(c) stem from simulations on lattices with 8x8 unit cells 
(L=8). In Fig. (a) both L — 6 and L = 8 results are included. 
As app arent and due to the inclusion of the magnetic field (see 
section III B ) size effects are next to absent until the magnetic 



length scale exceeds the lattice size (see text). 

Our results for the real part of the optical conductivity 
and associated resistivity p(T) = 1/(t((jj = 0, T) are plot- 
ted in Fig. ||. Here we consider J/t — 1.6. At high tem- 
peratures, T/t = 1/2, the data consists of a single Drude 
feature. As the temperature is lowered to T/t ~ 1/10, the 
weight under the Drude peak is reduced, and a feature 
at (jj/t ~ 1 is formed. A rise in the associated resistivity 
as a function of decreasing temperature is observed in 
this temperature range. This rise in resistivity originates 



from singular spin flip scattering on impurity spins and is 
reminiscent of the single impurity physics. By further re- 
ducing the temperature to T/t ~ 1/20 the Drude feature 
becomes sharper thus leading to a drop in the resistivity. 
Finally, at our lowest temperatures {T/t = 1/30,1/60) 
there is a marked upturn in the resistivity. At those 
temperatures the optical weight under the Drude peak 
starts to shift to frequencies cj/t ^ 0.1. 

J/t = 1.6 



X 



X 




Solid line: L = 6 
Dashed line: L = 4 



0.2 0.3 



0.5 



FIG. 9. (a) Uniform and staggered (inset) spin suscepti- 
bilities (b) charge susceptibilities and (c) specific heat for the 
DKLM at J/t = 1.6. 

The overall form of the resistivity versus temperature 
curve (Fig. ||a) allows us to define three temperature 
scales. The highest energy scale we will discuss corre- 
sponds to the resistivity minimum, Tc- For the consid- 
ered couphng Tc ^ 0.5t. We will define T* ~ O.lt as the 
temperature scale at which a maximum in the resistivity 
is observed. The magnetic scale T5 ~ 0.025i is defined by 
the energy scale below which the p{T) follows an activated 
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behavior. Next, we consider thermodynamic properties 
and concentrate on their behavior at the above defined 
energy scales. 

Fig. H plots the specific heat Cy, charge and spin uni- 
form susceptibilities (xc, Xs) as well as the staggered 
spin susceptibility {XsiQ)i Q = ^i/2) as a function of 
temperature at J/t = 1.6. From the technical point of 
view, the specific heat is computed using the recently 
proposed ME based method of Ref. |20j . As apparent in 
the temperature region Ts < T < T* one observes the 
following features, i) As indicated in Fig. ^ the specific 
heat is consistent with a 7T law. ii) the uniform spin 
and charge susceptibilities after going through a broad 
maximum seem to saturate. Both points are the charac- 
teristics of a Fermi liquid. The staggered spin suscepti- 
bility essentially measures the magnetic length scale at 
the considered wave vector. The inset of Fig. ^ shows a 
marked increase in this quantity at the spin scale T5. At 
the same energy scale, a sharp peak in the spin suscepti- 
bility is observed. Following the size effects, the data is 
consistent with a saturation of this quantity in the limit 
of zero temperature. At Tg the charge susceptibility ||b 
decreases rapidly. Upon analysis of size effects the data 
is consistent with the vanishing of this quantity at T = 0. 
Upon close analysis, an anomaly in the specific heat is 
apparent at T5. Thus, Ts marks the onset of substan- 
tial spin fluctuations which triggers the formation of an 
insulating state. That is the vanishing of the charge sus- 
ceptibility and an activated behavior of the resistivity 
below Ts- 



J/t = 0.8 



J/t = 0.8 



b 




FIG. 10. Optical conductivity as a function of temper- 
ature at J/t — 0.8. The inset plots the resistivity versus 
temperature. For those simulations the lattice size is L = 6. 




Solid line: L = 6 
Dashed line: L = 4 



0.1 0.2 , 0.3 0.4 0.5 

i /t 

FIG. 11. Same as Fig. | but for J/t = 0.8. 



Next, we consider smaller couplings, J/t — 0.8. Trans- 
port and thermodynamic properties are plotted in Figs. 
|l^and|ri|. As apparent from the conductivity data (Fig. 
no), no maximum as a function of temperature is ob- 
served. After hitting a minimum at Tc ~ 0.2t the re- 
sistivity grows monotonically as a function of decreas- 
ing temperature. The thermodynamic data of Fig. ^ 
equally shows no sign of the T* scale. That is, there is 
no energy scale for which the specific heat has a linear 
in T dependence. The only scales which are apparent in 
Figs.|l^ and |ll| are resistivity minimum Tc and the spin 
scale Ts- At Ts an anomaly in Cy is detected, a drop in 
Xc is observed and Xs exhibits a sharp peak. The peak 
in Xs is again located at the temperature scale at which 
XsiQ) shows a sharp increase. 

Our QMC results are summarized in Fig. |l2|. The 
highest energy scale, Tc, corresponds to the resistivity 
minimum, triggered by enhanced spin flip scattering off 
the impurity spins. As already observed for the KLM, 
this energy scale in a flrst approximation tracks J. At 
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high temperatures, it is tempting to try to understand 
this energy scale in terms of the single impurity Kondo 
model. In this case, we can use the closed form for the 
impurity resistivity (see and references theirin): 



R 



■imp{T) — 2 



HT/Tk) 



[\n\T/TK)+T^^S{S + l)\ 



1/2 



(20) 



where corresponds to the Kondo scale, ^ 
^g-w/J g^jj^fj ig i]^Q bandwidth. Assuming a power- 
law, T", for the resistivity of the conduction electrons, 
the resistivity minimum will scale as 



Tc - J^/" 



(21) 



at weak couplings. Given the above, our QMC result 
Tc ~ J follows by setting n — \ and we are left with 
searching for mechanisms which will lead to a linear in 
T resistivity. A possibility is scattering on spin fluctua- 
tions ll^lll. Note that both the KLM [|| and DKLM 
an enhancement of local antiferromagnetic spin-spin cor- 
relations between impurity spin and conduction electron 
((5g • S%)) is observed at an energy scale tracking Tc- 
Hence, magnetic fluctuations are potentially present at 
this energy scale. Clearly further work is required to jus- 
tify this hand waving argument, f] 




FIG. 12. Crossover scales for the DKLM as a function of 
J/t. Tc denotes the resistivity minimum. T* corresponds to a 
crossover scale at which a maximum in resistivity is observed. 
Below Ts at marked increase in the magnetic length scale is 
observed. We use the form e"^''""' to fit the T"*'"' and the 
form aT to fit the data for Tc- The crosses correspond to 
the Ko ndo t emperature of the single impurity problem. (See 
section 



[VC) 



Below Tc the resistivity grows as a function of decreas- 
ing temperature. For values of J/t > 0.8, the scale T* 



Note that in comparison to the KLM @, Tc is enhanced 
in the DKLM. Depleting further the DKLM of impurity spins, 
we should ultimately arrive to the single impurity case were 
Tc scales to infinity when the conduction band is modeled by 
free electrons. 



emerges. In particular and for the considered couplings, 
below T* the resistivity drops sharply, the specific shows 
a linear in T behavior, and both the spin and charge sus- 
ceptibilities seem to converge to finite values. In a first 
approximation, this is what expects for a Fermi liquid. 
Given the limited data, it is hard to pin down the func- 
tional form of T* as a function of J/t. Nevertheless, an 
exponential fit (e"^/""^) accounts reasonably well for the 
data. Due to the particle hole symmetry of the model 
which automatically leads to nesting, we cannot follow 
the metallic state to T = 0. At substantial spin-spin 
fluctuations set in and drive the system to an insulator 
with activated resistivity. For values oi J/t < 0.8, T* 
drops below Ts- Hence after a slow increase of resistivity 
below Tc, an activated behavior is expected below Ts- 



C. Comparing T* to the single impurity Kondo 
temperature 

It is now interesting to compare T* with the Kondo 
temperature of the single impurity problem: 

Hj = -tY, + + JS} ■ Sf.. (22) 

To obtain the Kondo temperature one can compute the 
impurity spin susceptibility: 



(23) 



where St{T) — e'^^' Sle - In the low temperature 

limit, T^ is the only scale hence Txi — F{T/T^) where 
is a universal function. 



Single impurity 




2.0, T/V* = 0.21 
1.6, T/V* = 0.12 
1.2,T,{/t = 0.06 



T/Tk 

FIG. 13. Impurity spin susceptibility for the single impu- 
rity problem. The calculations were done with the Hirsch-Fye 
impurity algorithm jl^ . For the conduction electrons we con- 
sidered a 512 X 512 lattice. 

In Fig. |l^ and for the three considered values of J/t 
a good collapse of the data is obtained with the listed 
values of the Kondo temperature T^. Comparison with 
Fig. ^ shows that the mean-field approximation produces 
next to identical results for the Kondo temperature. Re- 
call that the Kondo temperature for the impurity and 
lattice cases are identical in the mean-field approach. 

The important point however, is that within our ac- 
curacy T* compares very well with T^. This stands in 
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rough agreement with experiments on Ce^: Lai-^ Cug |Q. 
Here the single impurity Kondo temperature {x << 1) is 
given by ~ 3K. As x grows (0.73 < a; < 1), a maxi- 
mum in the resistivity versus temperature curve develops 
at energy scales in the range T"^ ~ 5 — 15K. It is also 
worth pointing out that below T* the spin susceptibility 
in CeCue seems to saturate Q. This stands in rough 
agreement with our results. Hence, T* corresponds es- 
sentially to the single impurity Kondo temperature 
below which the magnetic impurities are screened. It 
marks the onset on the heavy-electron state. 

V. CONCLUSIONS 

To summarize, we have presented a detailed numerical 
study of a depleted Kondo lattice. The depletion con- 
sidered here is very special, and corresponds, to the best 
of our knowledge, to the minimal regular depiction re- 
quired to obtain a two-dimensional metallic state in the 
J/t — > oo limit. 

From the technical point of view, the QMC approach 
is flexible and can be applied to any depletion pattern, 
random or regular. The constraint to avoid the minus 
sign problem is the particle-hole symmetry of the con- 
duction band. Hence various subjects may be studied. 
As function of depletion concentration, a metal-insulator 
transition should occur as observed in Ce3Bi4Pt3 when 
Ce is replaced by La In the strong coupling limit 
one expects this metal-insulator transition to occur at 
the percolation threshold. Furthermore, in the insulating 
phase random depletion of the /-sites should trigger mag- 
netic ordering between the magnetic moments induced by 
the depletion. Hence order by disorder is equally a phe- 
nomena which can studied in this model. To control size 
effects, we have introduced a magnetic field perpendicu- 
lar to the two dimensional lattice. We have shown that 
by scaling the magnetic field as BL^/^q extremely good 
size scaling to the thermodynamic limit is obtained. In 
particular already on small lattices the temperature scale 
at which size effects set in drops by an order of magni- 
tude when the magnetic field is included. This is the key 
to our results since it opens the energy window at which 
coherence effects may be studied. 

The QMC results for the DKLM show three temper- 
ature scales which are summarized in Fig. ^ i) The 
largest scale, Tc, corresponds to the minimum in the 
resistivity and tracks J/t. ii) Ts corresponds to the en- 
ergy scale below which the magnetic length scales grows 
presumably exponentially with decreasing temperature. 
Due to particle-hole symmetry and associated nesting 
properties those spin fluctuations drive the system to an 
insulating state. The ground state of the model remains 
a magnetically ordered insulator irrespective of the value 
of J/t. As in the KLM Q our data supports Ts - 
at weak couplings, in accordance with the RKKY inter- 
action. In the strong coupling limit Ts is exponentially 



suppressed and vanishes at J/t = oo. Note that in the 
absence of nesting one would expect the occurrence of a 
quantum phase transition between magnetically ordered 
and paramagnetic metallic states at a finite values of J/t. 
This is unfortunately not accessible to the QMC method 
due to severe sign problems, iii) The temperature scale 
T* as defined by the maximum in the resistivity versus 
temperature cure becomes apparent when it exceeds the 
spin scale Tg. This happens at J/t > 0.8. In the tem- 
perature range T* < T < Ts we see signatures of Fermi 
liquid character: the specific heat is roughly linear in T 
and the spin and charge susceptibilities appear to con- 
verge to finite values. Unfortunately, we cannot follow 
this metallic phase down to T = because of the onset 
of spin fluctuations below Ts. iv) We have computed the 
corresponding single impurity Kondo temperature (T^) 
with the Hirsch-Fye QMC algorithm. Within the consid- 
ered small intermediate coupling range, our results sug- 
gest ~ T*. 

Comparison with the mean-field results (see Fig. ^) 
shows that the spin scale is well reproduced. Clearly, the 
Tc scale is absent at the mean-field level. The mean-field 
Kondo temperature Tx is independent on the impurity 
density and hence is identical to that of the single impu- 
rity model. Comparison with exact QMC results shows 
that this mean-field approximation yields next to exact 
values of the single impurity Kondo temperature. Hence 
the mean-field Kondo temperature which marks a phase 
transition in this approximation evolves to the crossover 
scale when the model is solved exactly. 

Finally we comment on the coherence scale. At first 
glance it is tempting to associate the coherence scale with 
T* since signs of Fermi liquid character start emerging at 
this energy scale. However this identification will lead to 
contradictions. In particular, in the strong coupling limit 
the coherence scale scales to a finite value f\ but T* is 
expected to track T^ oc J. Hence we can only speculate 
that the coherence scale - as in the mean-field approach 
- lies well below T* and is masked by T5. 

We thank HLRS (Stuttgart) for generous allocation of 
CPU time on the Cray-T3E. I acknowledge useful discus- 
sions with S. Capponi, M. Feldbacher and A. Muramatsu. 
The Deutsche Forschungsgemeinschaft (DFG) is thanked 
for financial support under the grant number AS 120/1-1. 



^ In the strong coupling limit, the model reduces to the three 
band free electron model with Fermi line plotted in Fig. 
The coherence temperature is then given by the coherence 
temperature of the free electron model which is finite. 
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